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Abstract 

A system of interacting qubits can be viewed as a non-i.i.d quantum information source. 
A possible model of such a source is provided by a quantum spin system, in which spin-1/2 
particles located at sites of a lattice interact with each other. We establish the limit for the 
compression of information from such a source and show that asymptotically it is given 
by the von Neumann entropy rate. Our result can be viewed as a quantum analogue 
of Shannon's noiseless coding theorem for a class of non - i.i.d. quantum information 
sources. 
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1 Introduction 



In this paper we study the issue of compression of information from a particular class of quantum 
information sources, formed by systems of interacting qubits [see Section ^ for details]. Our 
aim is to quantify the minimal physical resources necessary to store the output from such a 
source or to transmit it through a noiseless channel. We shall use the words message, signal 
and output from a source interchangeably. The parameter that we minimise is the dimension 
of the Hilbert space to which a typical signal can be projected (i.e., "compressed") with high 
fidelity. In addition, it is expected that the interaction between qubits in the systems under 
consideration yields highly-entangled states; this is a motivation for the present work, even 
though the issue of entanglement is not discussed here. 

The analysis that follows shows that the data compression limit for output from such a 
source is given by the von Neumann entropy rate. This result can be viewed as a quantum 
analogue of Shannon's noiseless coding theorem for our class of non-i.i.d quantum sources. 
It can be considered as an extension of Schumacher's coding theorem '0 . 

Shannon's noiseless coding theorem quantifies the extent to which one can compress the 
information being produced by a classical information source. A standard model of such a 
source is described by a sequence of random variables Xi, X2, . . . , whose values Xi, X2, . . . , x„ 
represent the output of the source. For simplicity, consider random variables which take values 
from a finite alphabet of symbols or letters (extensions to infinite alphabets also hold). Let 
X_ := (Xi, X2, . . . X„) denote the sequence of random variables representing the source and 
x:= (xi, a;2, . . . Xn) the values that it takes. The source is described by a set of probabilities 

p(x) := Prob (X = x). 

An i.i.d classical source is one for which the random variables Xi, X2 . . . X„ are independent 
and identically distributed. In this case 

p(x) := p{Xi)p{x2) ■ ■ ■ p{Xn), 

where {p{x)} is the single symbol distribution. 

The main ingredient of Shannon's noiseless coding theorem is the Shannon entropy given 

by 

H{X) := -5Ip(^)log2 P{xy, 
for an i.i.d. source this reduces to H{X) = nH(X), where 

X 

In classical information theory one encodes the signal from a source into a string of binary 
digits (or bits). For purposes of storage and transmission, the aim is to encode the messages 
with sequences that are as short as possible. 
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An information source (classical or quantum) has redundancy, in the sense that certain 
outputs occur more frequently than the rest. This fact can be used to compress the source 
output: data compression is achieved by assigning shorter descriptions to the most frequent 
outputs of the source. The compression of data from a classical information source works as 
follows^ : A compression map, C", of rate R takes a sequence x = (xi, . . . , x„) of length n to 
a binary string of length [nR] (the symbol [■] denoting the integer part). A decompression map, 
D"', takes a binary string of length [nR] to a string of symbols of length n. The compression 
scheme is said to be reliable if with probability approaching one, as n — >■ cxo, {C^{x)) = x. 

Shannon's noiseless coding theorem indicates how well such a compression scheme works. 
More precisely, it asserts that for a large class of sources (i.e., stationary and ergodic), the mean 
length of encoded bit sequences is asymptotically given by the Shannon entropy, H{2C). More 
precisely, the data compression limit, which is the limiting number of bits per symbol, is given 
by the Shannon entropy rate: 

h:= lim - i7(Xi . . . X„). 

n— >oo fi 

An attempt to represent the source using fewer bits than this would result in a high probability 
of error when the information is decompressed. Hence a compression scheme of rate R is reliable 
only ii R> h. 

A quantum information source is defined in this paper by a set of distinguishable quantum- 
mechanical states \ipj)-, i.e., orthonormal vectors from a given Hilbert space, and a set of corre- 
sponding probabilities {k^}. We interpret the |'?Aj)'s as signals of the source and the k^-'s as the 
probabilities with which the signals are produced. Such a definition arises naturally from the 
density matrix formalism where a quantum-mechanical system is described by a convex linear 
combination of pure states: 

Here the |'?/'j)'s are identified as the orthonormal eigenvectors of p. The eigenvalue of p corre- 
sponding to I'j/'j) is Kj, and we have 

Kj > and ^ Kj = 1.0 

j 

More precisely, we deal with a sequence of 2" x 2" density matrices pn, n oo, and relate 
asymptotic properties of their eigenvalues k^J^^ to the von Neumann entropy rate. The von 

^More generally, one can consider a decomposition of p in terms of non-orthogonal states: 

P = ^Pi\(l)i){<l>i\, 

i 

where pi > 0, pi — 1. Such a decomposition is non-unique. Moreover, it is known that non-orthogonal states 
cannot be reliably distinguished. Hence, in the sequel, we focus on orthogonal decompositions. An exception 
is Theorem ^, which gives the limiting fidelity of the data compression scheme for general non-orthogonal 
decompositions . 
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Neumann entropy of a density matrix p„ is given by 

S{pn) = -trp„ log2 Pn = 4"^ log2 4"''' (1) 

and the von Neumann entropy rate by 

h= lim - S{pn). 

n—>oo 

A useful example is an i.i.d. case where p„ acts on a tensor product Hilbert space Tin = /C®"^ 
and is given by 

Pn = TT^". 

Here /C is a fixed Hilbert space (representing an "elementary" quantum subsystem) and vr is a 
density matrix acting on /C: 

i 

The eigenvectors iV'j"^) of Pn are tensor products 

\^f^) = \4>n)®\H)---®\^jJ- 
and its eigenvalues are given by 

(n) _ 

tvj — qj-^ . . . qj^ . 

This provides a convenient identification of label j as a "classical string" (ji, . . . ,jn) which will 
be emphasized by the notation j below. The von Neumann entropy is in this case S{7r^^) = 

In the case where dim /C = 2, the space 7i„ represents a system of n qubits. In analogy 
with classical data compression, it is desirable to represent typical outputs, l^/']"^), by vectors 
from a lower dimensional Hilbert space, thereby reducing the number of qubits needed for the 
source description. 

In his seminal paperS, Schumacher proved that the number of qubits necessary to represent, 
reliably, the signal from an i.i.d quantum information source is asymptotically given by the von 
Neumann entropy. More precisely, there exists a reliable compression scheme of rate R only 
when R > Sin) (under a suitable definition of fidelity). Schumacher's approach was developed 
further in o Ell. Extensions of Schumacher's theorem to some classes of quantum sources with 
memory have been established by Petz et al (see Bill and references therein). 

As was said before, in this paper we consider data compression for a class of quantum 
information sources which are modelled by a system of interacting quantum spins. This is 
an example of a quantum system with a strong coupling between the spins and with the 
environment and it does not fall into the classes of sources considered in the literature before. 
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Besides, we consider properties of eigenvalues which hold asymptotically with probability 
one; this is a refinement of results obtained inEJEJ. From the probabilistic point of view, our 
result is an analogue of the Shannon-McMillan-Breiman theorem (which is a version of the Law 
of large numbers), seefl. 

Even though we consider so-called quantum spin systems as models of a quantum source 
in this paper, our results also hold for sources modelled by quantum lattice gases, where the 
statistics (Rose or Fermi) of the particles is taken into account. 

Models of quantum information sources, based on large systems of interacting spins or 
particles, are being used increasingly in experiments with entanglement i'lll, as well as in 
theoretical research 111 . As mentioned before, our main result can be viewed as an extension 
of Schumacher's coding theorem to this class of sources. Section ^ contains a mathematical 
description of the class of systems under consideration. In Section |^ we prove that the data 
compression limit for such a class is given by the von Neumann entropy rate [see (|ll])]. The 
proof of the main theorem, which yields the data compression limit, is given in Section ^. 



2 Quantum spin systems 

We consider a quantum-mechanical system on a (i-dimensional lattice Z'^, with a spin-1/2 
particle attached to each site of the lattice. The particle can be either in an up-spin state 
(denoted by | |) ) or a down-spin state (denoted by | |) ). Hence, to each lattice site x E Z"^ is 
associated a Hilbert space JCx which is isomorphic to /C = C^, the single-qubit Hilbert space. 
For any finite subset X G Z'^, the corresponding Hilbert space is given by 

Here, and below, \B\ stands for the number of elements in a finite set B. Furthermore, we 
denote by Ax the algebra of 2^^^ x 2^^^ matrices acting in Hx ^ the local observable algebra. 
To each site x of the lattice, we associate a variable jx G {1, —1} such that jx = 1(— 1) when the 
spin at a; is t (i)- A configuration uj\ in a finite volume A C Z'^ is an assignment {jx, a; G A} of 
jx to each x E A; the set of configurations {uja} provides labels for a quasiclassical basis {IciJa)} 
in TiIa = (C2)®l^l. 

The physics of the system is described by an interaction, $ = which is a map taking 

finite subsets X C Z"^ to (self-adjoint) operators $x from Ax] see H. We study quantum 
systems that are small perturbations of classical ones. That is, we consider interactions of the 
form $ = $0 + Q, with 

= %x + Qx, (2) 

where, for all X, $ox is diagonal in the quasiclassical basis {Icja)} and Qx is small in norm 
(see below). We will write $o = {*&ox} and Q = {Qx}- 

The corresponding Hamiltonian H\ = J2xcA of a system confined to a finite volume 
A C Z'^ is written as a sum 

i^A = HoA + Va, (3) 
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where Hqa '■= J2xca^ox, Va '■= J2xcaQx- We make the following assumptions: 

(i) We consider translation-invariant interactions (for details, see®) i.e., $x+a — ^x, for 
all finite X C Z'^ and a e Z'^. The range of the interaction is defined as the supremum of the 
diameters of sets X from {X G Z'^ : X 3 and $x 7^ 0}. We use the £°°-diameter 

diam M : = max max \xi — vA , (4) 
and consider $ to be of a finite range, i.e., with R < oo. 

(a) The classical part $ox of ^x can be considered as a real- valued function on the set of 
configurations Ux in X ( i.e., an assignment {jx,x G X}). It is convenient to think of ^qx 
as a function of the infinite-volume configuration w = {jx,x G Z"^}, which depends on its 
restriction wx only. Similarly, Hq\ is a real- valued function of u depending on ua only. We call 
an infinite- volume configuration u = {jx, x G Z'^} periodic if jx = jx+a(^)y i = 1, . . . , rf, for all 
X E Z^ and a given collection of periods a^*-* = n^*-* e'-*^ where e*^*-* = (0, . . . , 1, . . . , 0) G Z*^ (entry 
1 at position i) and n^^^ is a given integer. A periodic co is called a ground state configuration 
for $0 if 

liminf ^ [iJoA(^') - Hoa{uj)] > 00, 
A/Zd |A| '-' 

for any infinite-volume configuration ut'. 

We assume that $o has a /inzte number of periodic classical ground states, a^^\ . . . ,(T^"^\ 
and satisfies the so-called Peierls condition 0. The latter is a condition for stability of the 
ground states relative to "local" perturbations. (See Hi and references therein for details.) 

(Hi) The term V\ is a quantum perturbation (i.e., [Hq/^, Va] ^ 0), with 

\\Qx\\<c\'^''^ 

for some constant c and some < A < 1. Here siX) denotes the number of sites in the smallest 
connected subset of the lattice containing X. We consider A as the perturbation parameter. 

Assumptions (i) - (Hi) constitute the framework of the so-called quantum Pirogov-Sinai 
theory [|, |, H, 0, 

We fix a boundary condition outside volume A, i.e., assume that the configuration on A*^ : = 
Z'' \ A coincides with a fixed reference configuration a, which is one of the periodic ground 
states aW,...,cr('") of i/oA- 

Since the interaction is of a finite range, the spins in A interact only with those spins in A^ 
that are in the envelopping volume A^: 

A'^ = {i G A'^ : dist(i, j) < R for somej G A}. 

^ Here and below, the symbol A Z'' is used for the thermodynamical limit, taken along a sequence of 
growing finite volumes Ai C A2 C . . . C Z'* of "nice" shape (e.g., hypercubes A„ = [— n, n]'^ n Z''), with 
U„A„ = Z'^. 
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Let be the orthogonal projection onto the subspace C T-^auas of dimension 2'^', spanned 
by states for which the configuration on is fixed to a. Then the Hamiltonian governing the 
spin system in A under the boundary condition a is given by 

X: 

The spin system with Hamiltonian if^ can be viewed as a system of interacting spins entangled 
with its environment. It is considered at a finite but low temperature. Due to the interaction 
between spins, the density matrix cannot be written as a tensor product of the density matrices 
of the individual spins and hence the quantum information source is non-i.i.d. The density 
matrix is written in the standard Gibbsian form: 

where /? > is the inverse temperature. The denominator on the RHS of is the partition 
function: 

The expectation of an observable A G ^auA' in the Gibbs state p'^'^ is given by 

(A)X = tr^. p-'^ A = ^tr«. A exp {-^HD . (6) 

Here and below, the trace is taken in the space 7i^; for notational simplicity, the subscript 
will often be omitted. For A := exp{iTH\), where r G IR, yields the characteristic function 
for the eigenvalues of the Hamiltonian H^: 

^-^(r) := {e-^^l. (7) 

The eigenvalues k^'^ of p°"'^ can be written as 



a,A 



(^f |e-^^X|V;f ) = _i_ exp {-Pi^pfm^pf)) , (8) 



where I'lpi''^), ■ ■ ■ , \'4'2w) ^-re the orthonormal eigenvectors of p"'^ (sometimes denoted by ipi'^ 
. . . , ip^w). The eigenvalues kJ'^ satisfy 



o-,A 

^ r 

,a,A 



E'^f = l- (9) 



The von Neumann entropy of p"' is given by 

^(p-'^) = -trp"'^ logs P"'^ 

-^4'"^ log2 4'^. (10) 
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The von Neumann entropy rate in this case is defined as 

A/Z-^ |A| 

= /Jet, (s'"' - /) (11) 

where Cq = log2 e and / and g^'^^ are standard thermodynamical functions (the free energy and 
the infinite volume energy per lattice site): 



We see that the von Neumann entropy rate h is well-defined if the above limits, ([T2|) and ([13|) , 
exist. The following theorem, proved in@, states that these limits do exist for the class of 
quantum spin systems under consideration. 

Proposition 1 Under the above assumptions, for /3 large and A small enough, the limits ( ff^ j 
and i^T^ ) exist. 

Remark: In this paper we deal with a sequence of density matrices p°^'^, A y TJ^ , not generated 
by a single state of a quasi-local algebra (see e.g. 'i). This puts us in a context different from 
that considered e.g. inttll. Hence we need Proposition |1] to guarantee the existence of the von 
Neumann entropy rate. 

In view of (^), the eigenvalues kJ'^, 1 < j < 21"^!, can be interpreted as the probabilities of 
the system being in the states IV'J''^). Let V"''^ be the corresponding probability distribution 
and consider a random variable K"''^ which takes a value kJ'^ with probability 

The data compression limit is related to asymptotical properties of random variables K'^''^ as 
A / T/. 

3 Data Compression Limit 

The main result of the paper is the following theorem. 

Theorem 1 Under the above assumptions, for (3 large and A small enough, for all 6 > 
where l(-) denotes an indicator function. 



8 



Note that 



where ]Ep^(-) denotes the expectation value with respect to the probabihty distribution 
Hence, 

hm lE?;^ = h, 

and Theorem |l| gives a Law of large numbers for random variables (— logg K"'^). 

The proof of Theorem |I] is given in Section ^. Here we discuss some of its consequences. 
The statement of the theorem can be alternatively expressed as follows: V 5 > 

lim P"'^ f2-l^l(^+'5) < K"^^ < 2-\^\^^-^A = 1. (16) 
A/Z'i \ - - J 

In other words, V e > and for A large enough, the eigenvalues kJ'^ of p"'^ satisfy 

2-mh+5) < ^j,A ^ 2-l^l(^-'5) (17) 

with probability > (1 — e). That is, the eigenstates IV'J'^) that correspond to eigenvalues kJ'^ 
satisfying ([T7|) are those which occur most frequently. We refer to them as typical states (or 
more precisely, 5-typical states). Let M."s^ be the subspace spanned by such states: 

A^^'"^ := spanjl^/^J'^) : (17) holds} (18) 

and lA^^'^l denote the dimension of this subspace. The following lemma establishes the growth 
rate of 

Lemma 1 For all 6 > 

ita^los.|A<fl = ft. (19) 

Proof: We follow a standard information-theoretical argument (see e.g. 0). From dll) it 
follows that the probability of a state being 5-typical is at least (1 — e) in the limit A y Z'^: 

(5) 

liminf ^kJ'"^ > 1-e. (20) 
j 

where the sum J2 is over those j's for which kJ''^ satisfies (|17D , i.e., G Ai'^'^. From (^) 

j 

(and the definition (|18D of the set A^^'^) we deduce that V e > 0, 

(<5) 

1 - e < E '^J'^ ^ 2-l^l('^-'5) IMs'^ (21) 
j 
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Also, from and (jl^) we have 



(S) 

2-\Mih+s) |_^-,A| < ^ ^j,A ^ ^_ ^22) 
From dll) and (H) it follows that 

Since this holds for all e > 0, we conclude that 

hm sup ^ log2 lA^^'^l <h + 6, hm inf log^ {M^A > h - 6. (23) 
Moreover, since 6 is arbitrary, 

lim sup ^ log2 lA^^'^^l = lim inf ^ loga lA^^'^^l- 

Hence limA/z<* jx\ ^^§2 l-^^'^l exists and is given by ([I9|). 

Lemma 2 Consider a quantum information source described by the density matrix p' 

(j.A \ ^ (T.A I /(T,A\ / /(T,A| 



,(T,A . 



Let 6e the von Neumann entropy rate [see If R > h then there exists a reliable 

compression scheme of rate R. 

Proof: Since there are at most 21^'^ 5-typical states (see Lemma |I|), one requires at most [\A\h] 
qubits to uniquely identify a 5-typical state. The data can be compressed as follows: 

Map each 5-typical state IV'J'^) to a quasiclassical state \x) = \xi) ® \x2) . . . \x[\A\h]), 
where x is a binary string of length [|A|/i]: 



X = (xi,X2, . . . ,a;[|A|h,]) e {0, 1} 



[\A\h] 



Clearly, this can be done in a one-to-one fashion, enabling us to recover any (5-typical state. 
In other words, the information contained in |A| interacting qubits is compressed into [\A\h] 
non-interacting qubits, which can be later decompressed unambiguously. In the limit A /" Z'^ 
this scheme succeeds with probability one. Hence, the data compression limit, for the class of 
non-i.i.d. quantum information sources considered in this paper, is given by the von Neumann 
entropy rate h. 

The following lemma shows that a compression scheme of rate R < h is not reliable. 
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Lemma 3 Let S\ be any set of eigenstates {|'?/'J' )} of p"'^ such that 

\Sk\ = 2[l^l«l, 

where R < h is fixed. Then for any e > and sufficiently large A 

E '^f < (24) 

je5A 

Proof: The LHS of (p^ gives the probabihty that an eigenstate of p"''^ belongs to the set iSa. 
We can write it as a sum of the probabihty that a state belonging to S\ is 5-typical and that 
it is atypical: 

(5) (5) ' 

E 4'^ = E 4'^ + E 4'^; (25) 

ie^A jG^A je^A 

here the second sum on the RHS of (^) is over the atypical states in S\. Choose 5 > such 
that R < h — S and < S < e/2. In the limit A y Z*^, the probability of atypical states is 
negligible. By ([T7| ) the total probability of atypical states can be made < e. There are atmost 
2\-\MRl (5-typical states in the set S\, each with an eigenvalue < 2~^^^^'^~^\ Hence, the first term 
on RHS of ( p5D is bounded by 

2-\A\{h-5) 2[\MR] < 2-1^1^/2, 

which goes to zero in the limit A y Z'^. 

We conclude this section with a theorem giving giving the data compression limit and the 
limiting fidelity of the compression scheme for general (not necessarily orthogonal) decomposi- 
tions of p"'^. 

Consider any representation of the density matrix p"'^: 

o",A _(T,A I , CT,A\ / , o-,A| 

i 

where G 7i\ are arbitrary vectors of unit norm [not necessarily orthogonal or even linearly 

independent) and p"'^ > 0, Y^iPi'^ = 1- To apply the above data compression scheme consider 
an orthogonal projection H : — > C, where C is a subspace of Ti^ such that the vectors 
|H0°^'^) are either collinear or orthogonal for different i (some of them may be 0). 

If such a projection exists then, necessarily, the vectors spanning C are eigenvectors of p'^'^ 
and each non-zero vector is collinear to one of these eigenvectors. If we take C to be 

the subspace J^f^, spanned by the 5-typical states IV'J''^) of p""'^, then to each non-zero vector 
\Il(f)i''^) we can assign a quasiclassical state associated with the eigenvector IV'J'^) collinear 
to IH^"^'"^) . Here x is a binary string of length < [log2 (dim C)] -|- 1. In this case, the compression 
scheme can be represented by the two maps given below: 

E : 10-'^) I^J'^) where I^J'^) G A^^'^ (26) 
C : IV'J'^) ^ k^-''^) where x^^^ G {0, 1}"; r < [hg^ (dim C)] + 1. (27) 
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We use the symbols E and C for the maps (|26| ) and (^) to denote encoding and compression. 
Note that map C is one-to-one. Hence, the quasiclassical state \x^^^) can be decompressed 
unambiguously to yield the 5-typical state |V'J'^)- However, map E is not necessarily one-to- 
one. Consequently, the original vector cannot be recovered with certainty from the state 
I'l/'J'^). Hence, we consider the following prescription for decoding the state l^/'J'^) (denoted by 
the map D): 

D : ^ 10-^), 

where satisfies the relation: 

{<Pk^\4^j'^) = max {(p'^'^l^f^). 
The fidelity of such a coding-decoding scheme can be defined as: 

i^A:=Epr^0r^ini0r^). (28) 

i 

The fidelity takes values between and 1 and equals to unity only when all the states |H0°^'^) 
are correctly decoded. In the following theorem we show that Fa tends to unity as A y TJ^ . 



Theorem 2 (i) Choose C to be the space of 6 -typical states of p'^'^: 

C = M'^s^ := span{|V^J'^) : 2-l^l('^-'^) > /tj'^ > 2-l^l(''+^)}, 

where the \ipi''^) 's are orthonormal eigenstates of p"''^ and kJ'"^ are their corresponding eigen- 
values. Let H be the orthoprojection Ti^ —>■ C. The fidelity F\ of the map U, given by 1^^), 
approaches one: 

hm Fa := }imXpT\<PT''\m''') = 1- 

(a) If, for some subspace T> C H.^, the orthoprojection H.- Ti^ has fidelity tending to one 
then ^ 

lim inf -—-logo (dimP) > /i, 

A/Z-^ |A| ^ ^ / - ' 

where h is the von Neumann entropy rate. 
Proof: To verify (i), write: 

i i j 

= E(v^f I Epr^i'/>r^)(</'r^ii(v^J'^ ^ c)\^f) 

j i 

= Y^i^J'^'lp-'^l^fW^'^eC) 
j 

= E'^f Idmlog'^f -^I<^)-1' (29) 



12 



by Theorem ^ Property (ii) is checked in a similar fashion. 

Remark. The argument in the proof of Theorem 2 does not depend on the nature of 
the density matrix p^'^ or space 7i\. In a somewhat different context, a statement similar to 
Theorem 2 was established in (see also the references therein). 



4 Proof of Theorem J 

In view of (|), eq. ([T^) is equivalent to 



E-f 1 (|co/5(^f l^^lV^f ) + log, - /.) I < 5) = 1. 

This fact, together with Proposition and eq. (|TTp reduces the assertion of Theorem to the 
following fact: V5 > 



(30) 



where f/'-""-' is defined through (ITSp. 

Eq. (^) is a Law of large numbers for the random variables (?/'J''^|i7^|'?/'J'^) (with respect 
to probability distributions V"'^). In terms of characteristic functions, (^) is equivalent to the 
following lemma: 



Lemma 4 For (3 large enough and X small enough, for any t G IR the following limit exists: 

hm ^(^)(t/|A|) = e^*^'^' (31) 

where is defined through ^ and g^'^^ by (^Wj). 
Proof : 

From (^ and (0) we have that 

V>^'\t/\A\) = (e^*^A/l^l)X = ^ ^ (32) 

Henceforth, we shall suppress the superscript a from the notation and 2°"'^. 
Expanding e**^A/|A| ^^le RHS of {^) we obtain 



n=0 \ 1^*- 



=: 1 + E^- (33) 

n=l 
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Let us first estimate the term Ti. 



A 



it 



I 



E tr^. ($xe-^^'^) I /S^ 

X: 

it ^ ^ 1 



A 



where 

Now 
and 



MM jeA MM jGA 



0.^:= E 



:= (tr$xe-^^-)/S^ 



< ll^xll, 

where 11 ■ 11 denotes the Hilbert-Schmidt norm of the interaction $x- Let 



Co := max | |- 

Due to the finite range of the interaction, we have that 

#{X9j|<l>x^0,jGZ'^} = (2 
for any site j G Z*^. Hence, 

mf)l\< E ^K$x)aI<co (2^^)'<oo. 

It is known that for (3 large enough and A small enough, the following limit 

A/'Z''- 

and defines the infinite volume Gibbs state i' H . 
Moreover, 

(0o)z-= E i^ii^^y^ E i^i'^xr vjGZ^. 

X90 1^ I X3j 1^ I 
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The last equality follows from the translational invariance of the interactions. 

Further, by using methods of B it can be shown that for (3 large and A small enough, the 
following bound holds: 



{^x)l-{<^xY\ < ||$x||c(s(X))r(dist(X,9A)) 



(41) 



here OK denotes the boundary of the volume A, s{X) is the number of sites in the smallest 
connected set of sites containing X, and the function r(r), r > 0, obeys 



|r(r)| < exp(— cir), 

where ci > is a constant depending on (3 and A. Note that V does not depend on A. 
Using our assumptions on $, one can prove the following Cesaro convergence: 

hm Ti = hm ^ Y.{^^)1 = ^^(0o)z- 



(42) 



Ayz 



A/'Z'' |A| 



(43) 



To prove (|43|) consider A to be a finite hypercubic volume [— n, ra]'^ fl Z'^ and define a subvolume 
A as follows: 

A := {z G A I dist(z, j) > log,/(A) Vj G dA}. (44) 
Here /(A) = 2n + 1 is the linear size of the volume A. In the limit A y Z'^, we have: 



We can write 



Now 



|A| 



(c) 



|A| 
\dA\ 



|A| 



1 

►0 



y") |A| - |A| 







(d) dist(A,9A) — ^ oo. 



it 



it 



it 



jeA ' ' jeA\A 



A/Zd A ^ A/Zd 

' ' iGA\A 



|A| 



supJ(e^^)X| 



jeA\A 



Hence, from (|38f) and 



(45) 



(46) 



(47) 



RHS of (47) < lim coltl (2^^) 

A/Zd ^ ^ 



d |A\ A| 
lAI 



0. 



(4J 
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Consequently, in the infinite volume limit, the second term on the RHS of (^61) goes to zero, 
thus allowing us to concentrate on the first term alone: 



it 



lim Ti = lim ^y{ef)i 



lim 



ill E H)! - (0o) 



A 



where 



jeA 
it 



jeA 



The last line of (|49|) follows from (^5| a). We shall prove that A = 0. Write 



A = lim 

A/Z'* 



jeA 



xcauaS 



jeA 



X3j 



:= + 

Recall that the interaction governing the system is of a finite range R. Define: 

Af ^ := {i G A|dist(i,j) < R}, J G A. 

Then we have 



Uil < lim 

A/"Z<* 



jeA xcA 



Using (|41|) we obtain 



'^^'-A>^.m^ E^^^II'^'xl|c(.(x))r(dist(x,9A)). 



jeA xcA 



(R) 



Set: 



We have that 



and for X C A 



C2 := sup c(s(X)). 



XcA 



(R) 



#{X|XcAf }=(2^^)', 
r(dist(X,9A)) < r(dist(A,(9A) - R). 



(49) 



(50) 



(51) 
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Hence, 

A/Zd |A| 



All < lim |t| 7nCoC2 i^'^^Y r(dist(A,M) -i?) = 0. (52) 



by (|3|d) and 

The second term on the RHS of (|50|) is bounded as follows: 



U2I < lim y J-|($^)-| 



< i>-.^2:E^I(*xn. (53) 

' ' jgA ' ' 

since A C A. Now 

#{X|X 9 j,X ^ A, J G A, <|.x ^ 0} = (22^)' 

and 

#{j G A| 3X 9 J, such thatX ^ A, $x ^ 0} = |Int(^)(A)|, 

where 

Int(^)(A) := {i G A|dist(i,M) < R} 
is the i?-interior of the volume A. We have 

It 



\A,\< hm^^ (2^^r co|Int(^)(A)|. 
k/z<i |A| ^ ^ 

However, 

|Int(^)(A)| < i?|9A|. 

Hence, 

\A,\ < limJt\\^Rco(2^^y 

A/Z<* |A| ^ ^ 

= 0, (54) 

by (H:). From (ID, (ID, (H) and (|5|) we readily get (||). 

This argument admits a generalisation for the n*^ term in the expansion on the RHS of 
(H). We have: 



n! \|A|_ 
1 / it' 



Y: tr($x,...$x„e-^^-)/H 



^A 



Xi ,...,Xn C 
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^' (|A|) Sn |X,r--|X„|^^^---^^"^ 



We prove below that for each n > 2, 



hm T„ = M!((0o)^.r. 
A/z-i n! 



Define volumes A*^") and A*^"); 



A(") = {{^l,...,^n)\^keA^l<k<n}, 

A(") = {(ji,...,j„)|(ji,...j„)gAW,L„>4")}, 



where 



Ln = L{ji, . . .,jn) ■■= min j^^min^^ [dist{jk, ji)] , min^ [dist(jfc, dA)]j 
and dA is the boundary of the volume A. The quantity is chosen so that 



hm £1"'' = oo 



and 



(a) 



|A(n)| 



|A{n) \A(n)| 



0, (C) ^TT^^ 0. 



I AW I ' I AH I ' ' ' I AW I 

[Note that |A(")| = \A\^.] Writing j = (ji, . . . , we prove ( p6D as follows: 



T 



1 

n! 



it 

W\. 



E(0t.---0t)A+ E (©^.---eDA 



jrGAt") 

T„(l)+T„(2). 



jeA(")\A 



(n) 



Now, 



QfXl < E ■■■ E 



1 



< 



(2 

< cS(2^'') 



■|(<l>x,...<fxJX| 
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Hence, 



|^{n) \ ^(n)| U 

lim TJ2) < lim ^ — -f^ — 



rid 



0, 



(62) 



by (|5P|c). Consequently, in the infinite volume limit, the only non-zero contribution to T„ arises 
from the term T„(l) on the RHS of (|60D . This term can in turn can be written as follows: 



1 I it 

lim TJl) < ^ \ m 
A/zd ^ ^ - n\ V A 



E [(0i---0l)X-((0o)z.r] + ((0o)z= 

jeA(") 



where 



it 



(63) 



B := lim T„(2) <-:,,, 



E •• ©Da -((eo)^.r 

■ • jGA(") 

The first term on the RHS of the last line of (|63D follows from (^a). We prove below that 



5 = 0. We can write B as follows: 

n 

I I 11 \ 

B : 



A/zd n\ (|A| 



E 

jeA(") 



E ■■■ E 

XlCAuA^ XnCAuA^ 



1 



1 



1^1 1 



\Xr, 



X 



((<i>x,...$xjx-($x....$x„r)+ E ■■■ E 4 



\x„ 



X (<l>xi . . . (^xXHXi (^AUA^ for some 1 < i < n) 
B, + B2. 



(64) 



By using methods of O it can be shown that for (3 large and A small enough, the following 
bound holds: 

71 n 
i=l i=l 

where 

A„ = A„(Xi, ...,Xn)= min{dist(Xi,Xj)|l <i< j <n}, 

and r(r) is a monotonically decreasing function of r, satisfying the bound (|4^). Further, recall 
that < Co and let C3 := sup^^^^ j^^^-gd c(s(Xi, . . . , X„)). For Xj C A and Xi 3 ji for 

1 < i < n, 

min {dist(Xi, X,)} > (4"^ - 2R). 

l<i<j<n 
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< 



nil* 



c(s(Xi,...,XOr(AO 



Hence, 

r(A„)<r(4")-2i?), 

and 



by (|8|) and (g^). Moreover, 



Itl" 1 

- ^ E ■■■ E ||$x,...$xJ|l(X,^AuA« forsomel<.<n) 

Itl" 1 

- }%\Wr\^ ^ E ■■■ E ||$x,...$xJ|l(X,^A forsomel<z<n) 
In fact, 

#{j G A(")| 3Xi 9 ji, . . .X„ 9 j„ such that Xi A for some 1 < i < n} = |Int(^)(A("))|. 
Here, as before 

Int(^)(AW) = {j G A(")| dist(j,9A(")) < R}. 
Hence, |Int(^)(A("))| < i?|9A(")| and 

|B,|< U„,^^e;fl(2«)"^ = 0, (67) 

' ' - A/zd |A|" n! ° V ; ' ^ ^ 

by (Hb). From (H), (|5l) and it follows that 5 = 0. Hence, from (H) and (H) one 
obtains 

hm T„ = M! ((0o)|,)" Vn > 1. (68) 
From (H) we now see, in view of Lebesgue's dominated convergence theorem, that 

oo (■j.\n 

hm y.(^)(t/|A|) = l + 5:^ ((0o)^.r. (69) 

The limiting energy density per lattice site, g^'^\ defined through (|13|), can be written as 

(7('^):= lim(ffA/|A|)X = lim V %^ = hm ^ V V (70) 
A^Z"^ ^' '^^ A/Z'^ |A| A/Z^|AI^ f^. |X| ^ ^ 
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Since the interaction $ = is assumed to be translationally invariant we can write 



Hence, (^) can be written as 
which proves Lemma ^. 
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